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ABSTR ACT 

A numerical treatment of acoustic waves in a Jet is described. The 
full time dependent Euler equations are used in both linear and nonlinear 
formulations. The computational region of integration is artifir ly 
bounded and boundary conditions are developed to simulate outgoing waves 
and to enable the computational domain to be substantially restricted. 
Higher order methods and coordinate transformations are introduced to 
further reduce the number of grid points as well as to increase the 
efficiency of the program. Numerical results are presented for time 
harmonic sources as will as for sources with more complicated time 
dependence. 
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I nt roduct ion 


This study is concerned with the numerical computation of time 
dependent acoustic waves in a Jet. The presence of inhomogeneit ies 
causes a significant change in acoustic behavior and aero-acoustic 
noise prediction schemes must take this into account in order to 
accurately predict aerodynamic noise. We will be concerned here only 
with the propagation of sound and not with its generation. Thus the 
acoustic sources will be modelled by arbitrary forcing terms in the 
equat ions. 

Previous workers (Schubert [1], Liu and Maestrello [2], Mungur et 
al. [3]) have approached the problem by deriving a convective wave 
equation for the pressure and then assuming tl it both tin; source and the 
solution have a time dependence of the form e*‘ jt . The use of this 
assumption in the convective wave equation leads to an elliptic equation 
for the acoustic pressure. In the case of zero mean flow this elliptic 
equation is the Helmholtz equation. 

A very efficient solution method for the elliptic equation was intro- 
duced by Schubert [1], He split the linear elliptic equation into two 
coupled nonlinear equations for the amplitude and phase of the time 
dependent solution. A motivation for this approach is that these quantities 
are smoother thar. the pressure and so can be efficiently approximated on 
a coarse grid. 

The meLhod used here differs from the above approach in two ways. 

First the time dependent equations are used rather than assuming a time 
harmonic behavior. This is because there are physical problems of interest 
where the source does not have a simple time harmonic dependence. For 
example, in a real jet the sources generally are convected downstream with 
the jet. Other sources of interest are those which have a pulse type time 



dependence. Numerical results for these types of sources will he presented 
here. 

A second difference from previous work Is that In the present study 
the first order equations of fluid dynamics are integrated. There are 
several limitations In the use of a convective wave equation for the 
pressure. The primary one Is that several physical approximations are 
necessary in order to eliminate the acoustic velocities (see [l]). These 
assumptions are generally of the type that the wave length of the sound Is 
small compared with the diameter of the Jet. Since only relatively low 
frequencies can he computed numerically the validity of these approximations 
is open to question. Furthermore, the exclusion of the velocities in the 
convective wave equation complicates the boundary conditions needed to 
simulate outgoing waves. Third the nature of the source term is less 
clearly defined than in the primitive equations of motion. For example, 
a source of mass or momentum can be easily included in the primitive 
equations but does not correspond to any simple type of source in the 
convective wave equation. On the other hand the primitive equations ran 
be easily modified to include nonlinear effects that are neglected in 
the standard linearizations. Finally explicit finite difference schemes, 
which are ideal for vector machines as the STAR-100, can be applied .o 
yield very efficient algorithms. 

The computational effort required to solve the first-order system is 
not substantially greater than that required for the convective wave 
equation. The second order wave equation requires at leaot two levels 
of storage while the first order system can be programmed with one level 
of storage for each of the unknowns. When the source has a harmonic 
time dependence then the tune dependent procedure can he considered as 
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an Iteration scheme to obtain the time harmonic solution. In this case 
further study is required to determine if this approach is competitive 
with the multi-grid methods used by Fix, f.unzburger, and Nlcolaldes [4]. 

The present method requires computational resolution of the solution 
on the basis of the wavelength. This will also be true of any method 
that solves the time Independent equations, either as a first order system 
or as a single second order equation. Such methods can not generate 
accurate solutions using grids as coarse as those used by Schubert. Indeed, 
the authors have found that a substantial fraction c* the total number of 
points us< 1 by Schubert is required merely to resolve the mean flow of the 
jet. 

The method of Schubert requires essent tally that the solution be composed 
of only one wave. Tf we define a wave as a function of the form 


w(r) = A(r)e 


ikS(r; 


where A and S are slowly varying real valued functions of the spatial 
point r, and k is the wave number, then A and S can be approximated 
with a coarser grid than the solution w. However, the assumption that the 
solution be composed of only one wave Is a very restrictive assumption 
which in general does not hold. For example, if the jet exits from a wall 
there will be reflections off t*e wall which will introduce other waves. Our 
results indicate that even for the problem considered by Schubert reflections 
occur because of the gradient of the mean flow. The method presented here 
permits the accurate computation of these cases but is restricted in the 
frequency range for which it can be applied. 
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I l . E quations 


The equations of mo' Ion for an inviscid axially symmetric flow in 
cylindrical coordinates are 

p t + (pu) z + * pv * r + ^ “ 0 

P z 

(2.1) u t + uu z + vu r + — - 0 

P r 

v +uv +vv + — -0. 
t z r p 

p is the density, u and v are the axial and radial components of 
velocity respectively and p is the pressure. We assume an Isentropic 
flow so that 

(2.2) p - Ap Y 

where A is a constant. These equations are to be Integrated for t > 0 

and for all space r > 0, -» < z < °°. In addition, we need to introduce 

additional terms Into (2.1) to simulate the acoustic sources. 

In order to incorporate properties of the mean flow we expand the 

T 

solution about the mean flow of the jet (p^.u^.v^ . 


(2.3) 



where t is a perturbation parameter determining the amplitude of the source. 
The acoustic approximation generally assumes e is small. In addition we 
non-dimensionalize the system by choosing new variables 
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where d is the diameter of the Jet, a is the ambient sound speed and 
p , p are the ambient pressure and density respectively. For simplicity 

■ eo ao 

we will drop the bars in the sequeal; all quantities will be non-dimen- 
sional ized. 

The mean flow for the Jet was obtained from experimental data 
(Maestrello [5]) which indicates that p^ is nearly constant and is, 

in the extreme case, only about 2 percent of u 0> Computational 
experience confirms that the Inclusion of v () makes little difference. 
Hence, for the remainder of this study we neglect v f) and assume that 
p^ is constant (p, B ). The mean velocity u 0 corresponds to a spreading 
jet with a spread of about 12°. For fixed z the mean flow has a 
maximum at the axis and decays with distance normal to the axis. Further 
details of this flow can be found in [5]. 

In order to facilitate the comparison of linear and nonlinear models 
we introduce new variables 

P - p' 

u * ( 1 + £ p ) u ' 

(2.4) 

v - (1 + e p )v' 

P - (( 1 + e p ) 7 -l)/e. 

For non-zero c the quantities u and v are part of the terms making 
up the non-dimensionalized momenta. The choice of dependent variables in 
(2.4) puts the resultant system of equations in conservation form. Upon 
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including a source term In the mass equation of (2.1). We can rewrite 
(2.1) as 


t)v 0 + v 

p t ♦ <pu 0 + u) z + (pv Q + v) f ♦ • F(z,r,t) 


(2.5) “ t + (ii(u 0 -H;u , )] z + |i(v 0 +iv')] r + j ( -^~ 
V + [v(u n +cu')l + (v(v 0 +Lv')l + i" 

L Y . 


0, r 0, r 


uv„ _ - vu„ _ + e u(F 4 ~ ) 


■ vv_ - uv_ + rv (F+ — ) 
0,z 0, z ' r 


0 < r < 00 , 


< z < ®, t > 0. 


If c is set to zero in (2.5) we obtain a linear system for the 
acoustic density and velocities. (The acoustic pressure is proportional 
to the acoustic density by linearizing (2.2)). For most of the study we 
have used the linearized equations, f * 0. Differences between the 
linear and nonlinear equations are discussed in section 5. 

Generally, we take the source as being of the form 
F(z, r , t ) - f (t)6(/ r 2 +( z-Zq) 2 ) for some axial point z^. The delta 
function source when implemer'.ed in the primitive equations of motion 
corresponds to a simple source of mass or momenta depending on which 
equation it is in. Axial dipoles and quadrapoles can also be handled 
easily. This is one of the advantages of working with the Euler equations. 
The delta function is modeled by a Caussian centered at the source point. 

In previous work ([l], [2]) the source was modeled by excluding a small 
circular region from the domain of integration and then imposing suitable 
boundary conditions. This approach was originally used for this study 
but was discontinued as being inconvenient for optimized vector computation 
on the CDC STAR- 100. 
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Kquatlons (2.5) must be supplemented by Initial and boundary conditions. 
The initial conditions specified are p, u, v ■ 0 at t • 0 (i.e. the source 
is switched on from a state of rest). Analytically the domain of integration 
is the entire plane; computationally we require a bounded domain with 
appropriate boundary conditions. This is treated in detail in section 4. 


III. Nu merical Scheme 

The system (2.0) can be written in the general form 


0.1) w t + F z + G r - H . 

We solve this system by splitting it into a sequence of one dimensional 
problems (see Strang [6], MacCormack [7]). Let \. , and denote 

the finite difference approximations to the one-dimensional equations 


w + F * H. 
t z 1 


w t + C r “ H 2 


respectively, where H - H ^ + Hj. Then, the solution to (3.1) is updated, 
with second order accuracy in time, by 


(3.2) 


w(t + 2At) * L L L L w(t). 

z r r z 


The solution was found to be insensitive to the decomposition of H 

v 

except that - had to be Included in L f . 
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For the numerical operators and L^, a spatially fourth order 

accurate extension of the MacCormack method developed by Gottlieb and 
Turkel [8) was used. The higher order accuracy was necessary in order 
to efficiently compute the solution at lar&u distances from the Jet 
exit. The higher order formulas allow a greatly reduced mesh without 
nnv loss of accuracy as compared with second order methods. Detailed 
comparisons of second and fourth order methods are given by Turkel [ 9 J . 

There it is shown that fourth order methods allow halving the mesh in 
each direction for a given error tolerance of about five percent. In 

addition the higher order method is about three times faster, to achieve 
comparable accuracy, as a second order scheme. In order to further reduce 
t iic number of mesh points an exponential coordinate stretching was intro- 
duced in both directions. This concentrates points near the axis and near 
the source. Finally at points where sufficient boundary conditions ar 
not prescribed a third order ext rapolat 'on of the fluxes is used as 
described in |9). 

The explicit nature of the fluxes F, G given by (2.5) allows their 
computation over the entire grid by vector operations. Hence, this 
algorithm is well-suited for a vector processor such as the STAR-100. 

The program was coded in 32 bit arithmetic in SL-1, a high level compiler 
written for the STAR. Including all factors, such as some recoding to 
facilitate vectorizat ion , the program ran 80 times faster on the STAR than 
on a CYBER-175 serial computer. 

IV. Bou ndary Conditions 

For computational purposes the infinite domain of integration for (2.5) 
must be replaced by a bounded domain on which appropriate boundary conditions 
must be specified. Since the infinite space problem is well posed for the 
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hyperbolic system (2.5), It Is necessary to develop radiation boundary 
conditions only to achieve accuracy and computat tonal efficiency. In 
principle, for any finite region and for a fixed time of integration one 
can place the artificial boundaries sufficiently far away so that the 
numerical solution Is accurate in the region of interest. This is not 
practical since in fact, one needs to constrict the region of Integration 
as much as possible, in order to reduce the number of grid points. Our 
experience indicates that this Is the fundamental computational difficulty 
in obtaining accurate solutions to (2.5). All this is fundamentally 
different from the Helmholtz equation where a radiation condition Is 
required to obtain a unique solution (cf. Hellwlg (10]). 

Two types of boundary conditi/rs i:re required corresponding to 
inflow and outflow situations. Inflow conditions occur where the jet 
enters the region of interest while outflow conditions occur at the 
far field cutoff where outgoing waves must be simulated. A complete 
description and justification of these conditions Is given by Bayllss 
and Turkel [11]. Here, a summary of the actual conditions, which were 
used in this study Is presented. 

A. Three types of inflow configurations are considered and are 
illustrated in figures 1-3. These figures are presented in cylindrical 
coordinates with axial symmetry understood. These Inflow conditions are: 

!• A semi-infinite pipe extending into free space (figure 1). 

2. A semi-infinite pipe extending from an acoustic baffle 
(figure 2). 

3. An anti-jet, that is a convergent upstream flow (figure 3). 
The anti-jet is obtained by reflecting the mean flow symmetrically about 
the plane of the nozzle exit. 
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Condition* 1 and 2 corre ond to physically realistic situations. 

In both cases two conditions must be Imposed at the pipe exit (dashed 
line in figures 1 and 2). We have Indicated (dotted line In figure 1) 
an outflow boundary directly above the nozzle; also included are far 
field boundaries (indicated by •) . Experimentally, one is not 

Interested In the acoustic field above the pipe even for case (1). 

The anti-jet is a nonphysical, computational device Introduced in 

the treatments of the convective wave equation in order to remove the 
necessity of specifying boundary condltons along the pipe and at the 
nozzle exit. These conditions are difficult to impose for the convec- 
tive wave equation since the velocities are not independent 
variables. However, the anti-jet formulation has several disadvantages. 
Besides being nonphysical it requires the solution of the equations in a 
much larger region than necessary. Finally, at the inflow (left) boundary 
u 0 is positive over portions of the boundary. Hence, two boundary con- 
ditions must be specified in order to have a well posed problem. We 
have not succeeded in developing such an additional condition. For the 
computational results computed with the anti -jet this inflow boundary was 
sufficiently far away so that u f) was smell. Using onlv the radiation 
boundary condition (to be described) instability was delayed long enough 
to permit a periodic solution to be generated downstream. Clearly the use 
of the anti-jet is unsatisfactory and inefficient for time dependent 
equations. The anti-jet was used In this study for only a few test cases 
for comparison purposes. All other cases used either boundary conditions 
1 or 2. 

For the inflow pipe we assume a constant mean flow u () and e ■ 0. 

We further assume that the flow is one dimensional, i.e. independent of 
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r. From characteristic theory one finds that the appropriate boundary 
condition)* at the Jet exit are 

(4.1) u ♦ p - 0 

(4.2) v • 0 . 

In |11] we present a Justification for these conditions for the full 
two dimensional problem based on an expansion of modes travelling; down 
the pipe; they help explain the numerical results presented in the next 
section. Rogers (121 has shown that these conditions imply that the 
energy flux through the pipe is negative, i.e., no energy flows from the 
pipe into the computational domain and hence (2.5) is well posed in the 
sense of Ho '..*•• i. 

B. The outflow boundary conditions must approximate the condition 

.or outgoing waves. The mean flow decays with the distance from the 

nozzle and when u^ is zero (2.5) can be reduced to a wave equation for 

p. The condition for outgoing waves can then be expressed in terms of 

9 ■ tan 1 — and d ■ (r^z")* 5 by 
z 



(see ’ll] and also Friedlander [13]). This implies that 

(4.4) p t + p d " "d + 0( V 

d 

which upon using (2.5) with u^ ■ 0, E. ■ 0 is equivalent to 
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(4.5) 


* P. 


5 ♦ o<i) 

d 


whore u‘* denotes the radial velocity in polar coordinate*!. Although 
these equations are strictly valid only when u () ■ 0 exterior to some 
sphere, we have applied It to very constricted domains with substantial 
success, as Is described In the next section. Computational experiments 
show that the correction in (4.5) is necessary In order to control 
false reflections at the boundary. The condition (4.5) is part ot 
a family cf boundary conditions asymptotic to any order In ^ 

(see [11]). 

The boundary condition (4.5) is coupled with cubic extrapolation 
on the fluxes for u and v as described in [9]. This enables us to 

solve for the variables on the boundary. Gottlieb, Gunzburger, and 
Turkel [14] demonstrate that numerical stability can be Improved by use 
of characteristic variables for those quantltltes which are not 
analytically prescribed. For the computations in this studv it was 
never necessary t'> use characteristic variables. Instead, the radiation 
boundary condition (4.5) was used to calculate the pressure after the 
two equations for u and v had been integrated. 

The radiation boundary conditions presented here differ from those 
introduced by Engquist and Majda[15]. Ours are asymptotic in the distance 
of the cutoff while the expansions of Fngquist and Majda are asymptotic 
in the deviation of the wave from normal incidence. The conditions at 
the nozzle (4.1 - 4.2) are formally equivalent to the first order condition 
in 115], the justification for which is based on the constraint v ■ 0 
In the pipe. An analysis of the appropriate equations shows that th ?se 
conditions are rigorously asymptotic in the distance of the pipe for 
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frequencies below a fixed limit, which includes all frequencies for which 
numerical studies are feasible. This is discussed in detail in [ill. 


V. Num erical Re sults 

The flrHt computations were designed to test the boundary conditions 
described in the previous section. Computations were performed comparing 
the nozzle inflow conditions (4.1 - 4.2) with the anti-jet. Computations 
with different frequencies as well as computations with pulse type forcing 
functions all yielded virtually indistinguishable results for both sets 
of boundary conditions. Therefore, all the remaining results were pro- 
duced for the configuration shown in figure 1. Here the outflow radiation 
conditions are imposed directly above the nozzle (dotted line in figure 1). 
This eliminates all upstream integrations. 

We first consider a time harmonic source. Then the forcing function 
in (2.5) is given by 


(5.1) 


F(z , r , t ) - H(t) 


■)[' 


cos wt 1 1 6(/ r" + (z-z () )" 


H ( t ) is a smooth approximation to the Heavisi .e function and z () is 2 
Jet diameters upstream of the nozzle. The linearized equations, f - 0, 
are solved. 

The quantity of physical interest is the relative sound pressure 
level (SPI.). Let p(z,r,t) be the non-dimensionalized pressure in (2.5). 
We then define 


(5.2) 


I (z.r) 


* 0 

2/ p (z , r , t )dt 
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Thin integral is taken over one period once a periodic oscillation has 
been established. If coswt in (5.1) Is replaced by e il and if we 
assume that 

p(z,r,t) - e iu)t p ( z , r ) 


then it follows that 


(5.3) I(z,r) • |p(z,r)|. 

We then define 

(5.4) SPL - 20 lo« 10 I (z,r) . 

Tills corresponds to the sound pressure level in the usual sense (see 
[ 1 ] ) . For general, nonharmonic problem the definition (5.2) is an exten- 
sion of the definition of (SFL) provided T is chosen large enough. 

The pressure is usually normalized relative to some reference pressure p () . 
Since only differences in SPL are important this normalization can be 
ignored. Differences between the SPL at different points are measured in 
decibels (d.b.). 

For the first set of experiments we consider the harmonic source (5.1) 
with different boundary conditions. The figures 4-7 display the SPL, 
relative to 90"\ at a fixed distance from the source and as a function of 
the angle 0. For these runs, unless otherwise stated, the exit Mach 
number Is .62, the non-dimensional frequency is 1.145, and the downstream 
boundary is at 65 diameters from the source. The radial cutoff is denoted 
by R^, and is measured in Jet diameters. 
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The fundamental physical effect to be expected here is the bending 
of sound waves away from the axis, a phenomena known as refraction. This 
results in a substantial reduction in the SPL as the axis is approached 
(i.e. 0 ► 0). This refraction effect decreases as 0 -* 90° so that the 
solution approaches the solution with no flow (i.e. the wave equation). 

This fact is used as a check on the code as well as a Justification for 
constricting the computational domain to a narrow region around the 
Jet axis. 

Figure 4a and 4b show plots of the SPL f or different values of R^,. 
When reference values at 90° are not available they are taken from the 
case with R^. ■ 57. It is clear that wc can severely constrict the domain 
normal to the Jet without seriously degrading the solution. The SPL is 
relatively Insensitive to polntwlse errors, however examination of the 
pointwiso solution indicates very small errors due to the constriction of 
the region. This is also true if the downstream cutoff is decreased. If 
the lower order term p/d is dropped in the radiation boundary condition 
(4.5) then much larger regions are required to control the reflections from 
the boundary. The solutions presented here agree closely with those 
obtained in [2] near the axis. The increase in the SPL that is found at 
mid-angles was found in all the runs with the harmonic source. This was 
not found in previous studies, probably because of a lack of resolution in 
the angular direction. 

In Figure 5 results are shown for the case of a Jet exiting from an 
acoustic baffle. In this problem there is a reflected wave from the 
baffle and it is clear from the figure that at higher angles there is 
at first a cancellation due to the interference of the waves and then 
a reinforcement. 
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In Figure b rorap.ir isons are presented with experimental results 
of Orande (wee |ll). We note that there is good qualitative agreement 
except that the rate of Increase of the SPL is steeper as 0 increases 
from get o. Comparisons with extensive experiments being conducted at 
NASA Langley Research Center will be presented subsequently. 

Computations were made with the nonlinear equations (2.5) with 
different values of l. Our results showed that the SPI. level was 
largely unaffected by the nonlinear terms except in the near field 
where the axial dip was slightly reduced for t .02 (based on a 

nondimensiona 1 P _ of 1). For higher values of c the scheme 
exhioited nonlinear instabilities which would require artificial 
dissipation to control. 

Hie time dependent calculations with the harmonic toning term 
Indicated that secondary reflections from the gradient of the mean flow 
are present near the axis. To demonstrate this more clearly we present 
numerical computations with a pulse type solution. In this case the 
source in (2.5) is given by 


(5.5) 


F(z,r,t) - f ( t ) a 


/ 2 "' 

. r + (z-z () ) 


where f(t) is a smooth approximation to a delta function. 

In Figure 7 the pressure is plotted as a function of time, at a 
fixed axial point 20 diameters downstream of the source. Results are 
presented for four different values of the exit Mach number. The striking 
feature to observe is the strong secondary wave occuring after the primary 
wave. We have verfied that this is Insensitive to numerical changes such 
as the position of the artificial boundaries and grid refinements. 
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The secondary wave is not a near field plienonmena and will occur 
at far field axial point* over a sufficiently long time Interval. The 
time of arrival of the secondary wave is inversely proportional to the 
exit Mach number while the amplitude is, to good approximation, propor- 
tional to the square of the Mach number. 

As a third example we consider time-harmonic sources converting down 
the axis of the Jet. We consider a model where sources drift downstream 
of the jet for one diameter and then are recreated at their initial posi- 
tion. This represents a reasonable model for the convection of acoustic 
* 

sources in a jet. Mathematically the forcing term (5.1) is modified to 
be 


(5.6) 


F ( z , r , t ) 


"(t)cosu) 


J~\ z-(z Q +a(t ))|‘ 


2 

r 


where 


(5.7) a (t) * c t (mod 1 ) . 

The convection speed c was chosen as .5 of the exit mach number of the 
jet and the frequency as 1.145. 

Analysis of the effect of the forcing term (5.6) on the solution of 
the system (2.5) is difficult even in the case of no mean flow where the 
acoustic density o satisfies the inhomogenous wave equation. Since the 
derivation of the wave equation from (2.5), involves taking the time deri- 
vatives of the first equation in (2.5), it is clear that the forcing term 
(5.6) will give rise to a combination of a dipole and a monopole as a 

This model was suggested by Lucio Maestrello. 
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forcing term for the wave equation. If the source convection speed c 
is reduced and the allowed length of drift in (5.7) increased, one expects 
a Doppler shift to a new frequency 

(5.8) w' ■ w[l + c*coe0) . 

This was in fact observed in the case of c very small and a large 
length of drift. 

We then consider more realistic parameters such as c ■ .3 and 
a drift length of 1. Figure 8 is a plot of the SPL change around circles 
centered at the origin of the source. The first thing to note is that 
these plots are very similar to those with no drift. There is a slight 
deepening of the dip as the axis is approached. This was observed in all 
tests with the drifting sources. Our results show that this deepening 
decreases as the convection speed of the sources increases and also that 
the field at angles close to 90° is basically unchanged. 

The fundamental difference due to the convection of the sources Is 
in the periodicity of the solution. Our results indicate that for the 
parameters of (5.7) the solution oscillates with the original frequency. 

The amplitude of each peak, however, varies substantially, and the solution 
exhibits a strong quasi-periodic behavior. 


VI. Conclusion 

The primitive (i.e. Euler) equations for the acoustic perturbations 
are solved for both t lie linearized and non-linear versions in axisvmmetric 
cylindrical coordinates. These equations are inherently more accurate, as 
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no approxlmat ions are made, as compared to the convective wave equation. 

There 1 h little computational disadvantage in solving the first order system 
as opposed to solving the time dependent convective wave equation. These 
equations are suitable for both harmonic time dependence and more complicated 
temporal behavior. 

Boundary conditions have been developed which permit the solution to 
be computed over severely restricted domains. These conditions are asymptotic 
In the distance of the artificial cutoff. The approximation (4.5) is 
easy to implement and allows very little reflection of the outgoing 
waves into the domain of Integration. The boundary conditions coupled 
with a fourth order numerical method permit the computation of accurate 
solutions with reasonable grids. The explicit code results in great 
efficiencies on a vector machine (see [16]). 

Solutions generated with a time harmonic source agree closely with 
previous computations in describing the dip in pressure near the axis. 

Differences away from the axis can probably be attributed to the increased 
resolution of the present computations. Further results are described 
with other time-dependent sources as pulse type forcing functions and 
moving sources. These reveal new phenomena which can not be calculated 
by previous codes. Extensive computations have indicated that the relative 

SPI. level is very stable under large perturbations. In order to dlstineulsh the 
effect of various sources one must examine the pressure as a function of 
time. Physical experiments are presently in progress to verify several 
of the computations presented in this study. 
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Figure 1. Computational domain for semi-infinite pipe. 
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Fipjre 2. Semi -infinite pipe extending from 
acoustic baffle. 
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Figure 4a. SPL (relative to 90°) for different values 

of R,p on a circle 51.4 diameters from source. 
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Figure 4b. SPL (relative to 90° for different values 

of Rrj. on a circle 19 diameters from source. 
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Figure 5. Comparison betwee jet exiting from a baffle 
and jet in free space (ft= 1.145, distance 51.4 
diameters from source). 
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Figure S. Comparison with experimental results of 

Grande (ft = 1.055, Mach number = .9, distance 
of 100 diameters from source). 
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Figure 7a. Pressure from pulse forcing term (Mach number ■ .01. 
distance of 20 diameters from source). 
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Figure 7b. Pressure from pulse forcing term (Mach number * .31. 
distance of 20 diameters from source). 
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Figure 7c. Pressure from pulse forcing term (Mach number - .62, 
distance of 20 diameters from source). 
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Figure 7d. Pressure from pulse forcing term (Mach number = .85. 
distance of 20 diameters from source). 
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Figure 8. Relative SPL level for sources 
converting down jet (T<= 1.145, 
Mach number = .62, convection 
speed - .31, drift length «|). 
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